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Maxwell’s equations cannot describe a homogeneous and isotropic universe with a uniformly 
distributed net charge, because the electromagnetic field tensor in such a universe must be vanishing 
everywhere. For a closed universe with a nonzero net charge, Maxwell’s equations always fail 
regardless of the spacetime symmetry and the charge distribution. The two paradoxes indicate that 
Maxwell’s equations need be modified to be applicable to the universe as a whole. We consider 
two types of modified Maxwell equations, both can address the paradoxes. One is the Proca-type 
equation which contains a photon mass term. This type of electromagnetic field equations can 
naturally arise from spontaneous symmetry breaking and the Higgs mechanism in quantum field 
theory, where photons acquire a mass by eating massless Goldstone bosons. However, photons 
loose their mass again when the symmetry is restored, and the paradoxes reappear. The other 
type of modified Maxwell equations, which are more attractive in our opinions, contain a term with 
the electromagnetic field potential vector coupled to the spacetime curvature tensor. This type 
of electromagnetic field equations do not introduce a new dimensional parameter and return to 
Maxwell’s equations in a fiat or Ricci-fiat spacetime. We show that the curvature-coupled term can 
naturally arise from the ambiguity in extending Maxwell’s equations from a flat spacetime to a curved 
spacetime through the “minimal substitution rule”. Some consequences of the modified Maxwell 
equations are investigated. The results show that for reasonable parameters the modification does 
not affect existing experiments and observations. However, we argue that, the field equations with 
a curvature-coupled term can be testable in astrophysical environments where the mass density is 
high or the gravity of electromagnetic radiations plays a dominant role in dynamics, e.g., the interior 
of neutron stars and the early universe. 

Keywords Classical general relativity. Maxwell’s equations. Classical fields in curved spacetime. 
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I. INTRODUCTION 

It has long been known that when some fundamental laws of physics that have been well established in labs are 
applied to the universe as a whole troubles or paradoxes may arise. A famous example that perfectly illustrates the 
above claim is the application of the second law in thermodynamics to the entire universe, which has led people to 
propose the idea of heat death for the ultimate fate of the universe. The heat death hypothesis states that, when time 
lasts long enough, the universe will reach a thermodynamic equilibrium (and hence maximum entropy) state with 
absolutely uniform temperature and energy, therefore any work process and ordered movement sustained by energy 
will not be possible any more. It was first proposed by Thomson (Lord Kelvin) and then popularized by Helmholtz 
and Rankine [3,0. However, after Thomson, Planck has put doubt about the heat death by criticizing the definition 
for the entropy of the universe and has pointed out that such a definition has no meaning (i ; see also Bi for 
comments on the concept of the entropy of the universe). 

Modern investigations on cosmology and gravity have put more doubts on the heat death. It is well known that 
structure formation and evolution in the universe is dominantly governed by the law of gravity. The heat capaci^ 
of a gravitationally bound system is negative, so that the temperature of the system increases as it looses energy [6|. 
Hence, it appears that the role of gravity is to keep the universe out of thermodynamic equilibrium 0. Einstein’s 
theory of general relativity has led people to believe that the ultimate fate of a massive gravitationally bound system 
is the formation of black holes. However, studies on quantum properties of black holes have revealed that black holes 
are thermodynamic objects and contain a huge amount of entropy pro portional to the area of their event horizons 
[1,0, and black holes evaporate by emitting thermal radiation [lol lllj[. For recent discussions on the entropy and 
structure formation in the universe based on the concept of gravitational entropy, please refer to [SO. Although 
the problem of heat death has not been ultimately addressed yet, the interplay of general relativity and quantum 
theory has refreshed our understanding about the concept of entropy and may provide a solution. 
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Another historical example for the inconsistency of a fundamental law of physics with a global universe is the 
so-called Neumann-Seeliger paradox [l^. In the 1800s, Neumann [T^ and Seeliger [T^ independently found that 
when Newton’s inverse-square law of gravity was applied to an infinite and static universe with a uniform distribution 
of matter, one would get that the gravitational field at any position in space had a magnitude of infinity. The 
inconsistency can be more easily seen by the following simple argument. For a universe with a homogeneous and 
isotropic distribution of matter, the gravitational field must be vanishing at every space position since otherwise the 
direction of the gravitational field g would specify a preferred direction at that position. If Newton’s law of gravity 
is written in the form of 


V • g = -47rGp , (1) 

where p is the mass density of matter and G is the gravitational constant, one immediately sees that it is violated 
since the left-hand side of the equation is zero but the right-hand side is nonzero. 

The same problem is encountered in astrophysics when people consider linear perturbations to a static and uniform 
system. For a gas of uniform density, pressure, and zero velocity, to the zeroth-order (i.e., for the unperturbed state) 
the equation of Newton’s gravity is violated. Astrophysicists steer by this problem by simply ignoring the unperturbed 
equation for Newton’s gravity and applying equation (ED only to the perturbed density and the perturbed gravitational 
field [i[i3. This treatment is usually called the Jeans swindle. 

The trouble in application of Newton’s law to the universe as a whole, or to a globally uniform system, is essentially 
caused by the fact that the entire Newton’s theory was based on the concept of absolute spacetime. In a universe 
that is homogeneous and isotropic, at every space position an object must suffer a zero external force so is an inertial 
frame. According to Newton’s first law, any two inertial frames must be at rest relative to each other, or move to or 
away from each other with a constant velocity. However, because of the existence of mass in the universe, any two 
inertial frames at a distance should accelerate to each other according to Newton’s law of gravity, unless a positive 
cosmological constant is introduced. This paradox constitutes a fundamental and conceptual problem for describing 
the dynamics of the universe with Newton’s theory.^ 

The difficulty of Newton’s law of g ravity in application to cosmology completely disappeared after Einstein invented 
the theory of general relativity [2a and Hubble discovered that the universe is expanding [2^. In fact, general 
relativity is the first and the unique widely accepted physical theory that can consistently describe the dynamics of 
the global universe. In general relativity, the usual concept of gravitational field is replaced by that of spacetime 
curvature: the presence of mass makes spacetime curved and hence affects the motion of objects in the spacetime. 
The equivalence between a local gravitational field and the acceleration of a point mass makes it possible to choose a 
frame at any position in a spacetime (the local inertial frame) in which the gravitational field vanishes. Einstein’s field 
equation can describe the dynamics of a homogeneous and isotropic universe without any problem, although it reduces 
to Newton’s equation in the limit of weak fields, slow velocity, and small scales [2^ [ 23 . So, the Neumann-Seeliger 
paradox does not exist in general relativity. The Jeans swindle can also be properly justified in the framework of an 
expanding universe, at least in principle (28|| . 

When Maxwell’s equations of electromagnetic fields are applied to the universe as a whole, a trouble similar to the 
Neumann-Seeliger paradox also arises. If we replace the g on the left-hand side of equation ED by the electric field 
E and the —Gp on the right-hand side by the electric charge density pe, equation ED becomes the Gauss’s law of 
electromagnetism, i.e., the first equation in the complete set of Maxwell’s equations 

V • E = 47rpe . (2) 

For a homogeneous and isotropic universe filled with a uniformly distributed charge, we have E = 0 everywhere since 
otherwise E would provide a preferred direction. Then Gauss’s law is violated, since the right-hand side of equation 
([2]) is nonzero. This indicates an inconsistency of Maxwell’s equations with a homogeneous and isotropic universe 
containing a net charge. One might object to the ratiocination by the arguments listed below, which we will show 
not to hold. 

1. The Gauss’s law in equation Q is an equation in a flat spacetime, but the universe is a curved spacetime. That 
is true. However, in Sec. mi we will show that in the framework of general relativity the above ratiocination also holds 


^ Although the dynamic equations for the evolution of and structure formation in a homogeneous and isotropic universe can formally be 
derived with Newton’s law of gravity and motion, it does not deny the fundamental and conceptual problems of Newton’s theory in 
application to the universe as a whole [T3. Il9l| . For example, due to the fact that in Newton’s theory light travels at an infinite speed 
and gravity propagates instantaneously, the formulation of the Newtonian equations applies only to a region smaller than the cosmic 
horizon. Extension to scales comparable to and larger than the cosmic horizon must include relativistic corrections sm. In fact, 
once the distance is extended to cosmological scales, the definition for distance in the Newtonian equations becomes ambiguous since 
there are multiple and distinct definitions of distance in an expanding universe [2^ . 
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and the inconsistency still exists. In fact, we will see that equation Q remains valid in the standard Maxwell theory 
extended to a curved spacetime, if the V is interpreted as a covariant derivative operator on a spatial slice of the 
curved spacetime. 

2. A universe must have a zero total charge then the inconsistency does not exist. However, there is no known 
first principle to force the total charge in a universe to be exactly zero. At least there is yet no such a principle that 
we know. In fact, a nonzero total charge can arise from several possibilities, including difference in the magnitude 
of electron and proton charges, nonvanishing charge of neutrons and neutrinos, and asymmetry of matter and anti¬ 
matter. Observational and experimental tests of the charge difference and asymmetry are being persistently carried 
on by researchers, although positive results have not been obtained yet |29l - l32l |. 

3. Maxwell’s equations determine that the total charge in a universe must be exactly zero. However, there is no 

reason to believe that Maxwell’s equations must be correct on cosmological scales. Although gauge invariance is often 
taken to argue for a massless electromagnetic field, in quantum field theory it is well known that a gauge field can 
acquire a mass through spontaneous symmetry breaking and the Higgs mechanism [H, [s^. So an electromagnetic 
field equation with a photon mass term is possible. Observational and experimental limit on the photon mass are 
being actively pursued by many researchers and references therein). 

4. The real universe is not precisely homogeneous and isotropic hence the above argument does not apply. If one 
considers linear perturbations to a homogeneous and isotropic universe with a uniform distribution of charge, similar 
to the case of linear perturbations to a static and uniform astrophysical system studied by Jeans one may get the 
correct first-order perturbed Maxwell equation by the Jeans swindle. However, the zeroth-order Maxwell equation 
remains unsolved. In addition, as will be shown latter in this paper, in a closed universe with a nonzero charge the 
Maxwell equation always fails regardless of the symmetry of the spacetime and the distribution of charge. 

The inconsistency problem of Maxwell’s equations in application to a homogeneous and isotropic universe was 
already encountered by Lyttleton & Bondi [37| in 1959, when they explored the physical consequences of a general 
excess of charge. They proposed that the observed expansion of the universe is caused by the repulsive force arising 
from a uniformly distributed net charge due to a tiny difference in the magnitude of electron and proton charges. 
They found that Maxwell’s equations must be modified to be consistent with the model of a steady de Sitter universe, 
and a steady creation of charge is required (see the comment by Hoyle [3^ and the response by Lyttleton & Bondi 
[^1. However, as Barry correctly pointed out, in a homogeneous and isotropic universe the electric field must 
be vanishing as forced by the symmetry of the spacetime, hence there cannot exist a repulsive electromagnetic force. 
More recent discussions on a charged universe and notes on the inconsistency with Maxwell’s equations can be found 
in Refs. [4ll - l4l| . and references therein. 

In this paper, the trouble and paradox arising from application of Maxwell’s equations to a homogeneous and 
isotropic universe with a net charge are presented and proved with a rigor analysis in the frame work of general 
relativity. They are also generalized to a more general case: in a closed universe or spacetime (i.e., the spatial section 
of the spacetime is a three-dimensional compact manifold) with a nonzero net charge. Maxwell’s equations always 
fail without any requirement on the symmetry of spacetime and of the charge distribution. Then, we investigate 
how to address the troubles and paradoxes by considering modifications to Maxwell’s equations. First, we show that 
an electromagnetic field equation with a photon mass term (i.e., the Proca equation) can address the paradoxes. 
Although spontaneous symmetry breaking and the Higgs mechanism can naturally generate a photon mass term in 
the electromagnetic field equation, we argue that it is not the ultimate solution since the photon mass disappears 
when the symmetry is restored. Maxwell’s equations with a native photon mass term is not favored because they 
introduce a new and unusually small dimensional parameter (a photon mass). 

Then, we consider to modify Maxwell’s equations in a more fundamental way by introducing to the field equation 
a term which couples the electromagnetic field potential vector to the Ricci curvature tensor of the spacetime. This 
is more attractive than introducing a photon mass term, since it does not introduce any new dimensional parameter. 
In a Ricci-flat spacetime, including the flat spacetime as a special case, the curvature-coupled term disappears and 
the electromagnetic field equation returns to the Maxwell equation. This means that in a Ricci-flat spacetime gauge 
invariance is restored. Since a non-empty universe must not be Ricci-flat, the paradoxes can be addressed by the 
inclusion of a curvature-coupled term in the electromagnetic field equation. If there is a cosmological constant in 
the universe, the curvature-coupled term in the electromagnetic field equation implies a photon mass determined by 
the cosmological constant. In this sense we can say that the new electromagnetic field equation implies the Proca 
equation. 

Besides the above mentioned inconsistency problem of the standard Maxwell equation in application to a universe 
as a whole, we have a stronger motivation for inclusion of a curvature-coupled term in the electromagnetic field 
equation. In general relativity, a law of physics in a flat spacetime (except the law of gravity, of course) is usually 
generalized to a curved spacetime by the “minimal substitution rule”. That is, in the law of physics, the Minkowski 
metric tensor is replaced by the general metric tensor in the curved spacetime, and the ordinary derivative operator 
is replaced by the covariant derivative operator associated with the general metric tensor . There is a well-known 
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ambiguity in this rule. If the law of physics contains second-order derivative operators acting on a vector or a tensor, 
the order of covariant derivative operators matters in a curved spacetime. In a curved spacetime, when a second-order 
derivative operator acts on a vector or a tensor, exchange of the order of derivative operators causes the presence 
of the Riemann curvature tensor in the equation. Because of this ambiguity in the “minimal substitution rule”, 
generalization of a law of physics in a curved spacetime is determined only up to addition of a term coupled to the 
spacetime curvature. Hence, it is natural to consider an electromagnetic field equation with a curvature-coupled term, 
since the field equation contains second-order derivatives of the potential vector. 

The “minimal substitution rule” is essentially reflection of the principle of equivalence, which states that at any 
point in a spacetime it is possible to choose a free-fall frame where in a sufficiently small region the law of physics 
takes the form as in a flat spacetime [13, • The critical point in the statement of the principle of equivalence is 

that the region in the question must be sufficiently small: it must have a size much smaller than the radius of the 
spacetime curvature, so that any term in the law of physics inherently coupled to the spacetime curvature can be 
ignored. Therefore, the presence of a term coupled to the spacetime curvature in the law of physics does not violate the 
principle of equivalence. Such a term cannot be recovered by the “minimal substitution rule”, since it does not show 
up in a flat spacetime. In a region that is not small compared to the curvature radius of the spacetime, the principle of 
equivalence breaks down and the effect of a curvature-coupled term can be important. For the case of electromagnetic 
fields, the presence of a curvature-coupled term in the field equation can represent interaction of electromagnetic fields 
with the spacetime curvature. In fact, when the standard Maxwell equation in a curved spacetime is expressed as 
a second-order differential equation of the potential vector (i.e., as a wave equation) a curvature-coupled term also 
appears 271, whose effects have been studied in [dsl. li^. 

In the paper we also discuss some consequences of the modified field equations, including the time delay in propa¬ 
gation of photons with different energy, effects on the statistical mechanics of photons, and constraint on the charge 
excess in the universe and the difference in the magnitude of electron and proton charges. We also investigate the 
testability of the new electromagnetic field equation with a curvature-coupled term. We find that, the effect of the 
curvature-coupled term can be detectable in an astrophysical environment with a high mass density, or in the early 
universe when radiation drives the dynamics of cosmic expansion. 

Throughout the paper, we use the geometrized and Planck units with G = c = h = 1, where c is the speed of light, 
and h is the reduced Planck constant. However, in a few places the units are restored to get the magnitude of physical 
quantities. 


II. FAILURE OF MAXWELL’S EQUATIONS ON COSMOLOGICAL SCALES 


As a physical law except the gravitational field equation in a flat spacetime is transplanted to a general curved 
spacetime, one usually adopts the “minimal substitution rule”: simply replacing the Minkowski metric rjab appearing 
in the law by the general metric gab, and the derivative operator da associated with rjab by the derivative operator Va 
associated with qgh. With this minimal substitution approach. Maxwell’s equations in a curved spacetime take the 
form (see, eg., ( 23 ]) 

= -dTT J'' , (3) 


and 


V[aF'bc] = 0 . 


(4) 


Here Fab is the anti-symmetric tensor of the electromagnetic field, J“ is the charge current density 4-vector, and 
square brackets in the index of a tensor denote anti-symmetrization of the tensor. 

In this section, we show that the inhomogeneous Maxwell equation (0) fails when it is applied to a homogeneous 
and isotropic universe with a uniformly distributed net charge. We call it Type I Paradox. In addition, if the spatial 
section of the universe is compact (i.e., the universe is closed), the inhomogeneous Maxwell equation ^ always fails 
if the universe has a nonzero net charge, without any requirement on the symmetry properties of the spacetime and of 
the charge distribution. We call it Type II Paradox. 

For a universe with a homogeneous and isotropic distribution of matter, which is considered to be consistent in 
high precision with all modern observations on our universe, we expect that the electromagnetic field must vanish 
everywhere. Otherwise, the electric field or the magnetic field will provide a preferred spatial direction at a position 
in the universe, which conflicts with the assumption that the universe is homogeneous and isotropic (i.e., the so-called 
cosmological principle 22|). Then we must have Fab = 0, and by equation ([3]), J“ =0 everywhere in the universe. 
This immediately implies that the inhomogeneous Maxwell equation cannot describe a homogeneous and isotropic 
universe with a uniformly distributed net charge. 
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FIG. 1: A compact spacelike hypersurface Et in a spacetime {M,gab) is divided into two parts of volume Vi and V 2 by a 
two-dimensional surface S. The surface S can be considered as the boundary of Vi with a normal vector n“. It can also be 
considered as the boundary of V 2 with a normal vector —n“. The electric field vector on the surface S is denoted as The 
total charge contained inside Vi is Qi. The total charge contained inside V 2 is Q 2 - Application of the Gauss’s law in equation 
([8)1 to the volume Vi and V 2 leads to the conclusion that Qi -f Q 2 = 0. 


The above statement can be more clearly seen if the Maxwell equation (|31) is converted to the form of Gauss’s law. 
To do so, consider an observer comoving with the expansion of the universe with a 4-velocity where t is 

the cosmic time. The electric field and the magnetic field measured by the observer are related to the anti-symmetric 
tensor F^b by 


Ea = FabU^ , (5) 

and 

Ba = -\eabcdU^F‘^'^ , (6) 

respectively, where Cabcd is the totally anti-symmetric tensor of the positively oriented volume element associated with 
the metric gab- Note, EaU°‘ = BaF^ = 0. 

On a spacelike hypersurface Et orthogonal to u°‘ (defined by t = const, such a hypersurface always exists for a 
homogeneous and isotropic universe), we can define a spatial metric on E( induced from gab by hab = gab + UaUb, and 
a spatial derivative operator Da associated with it. Then, from equation (|31) we can derive that 

DaE^^ = Attp, , ( 7 ) 

where pe = —UaJ°' is the charge density measured by the observer. For a universe with a uniformly distributed charge, 

Pe is a function of the cosmic time t. But, as discussed above, for a homogeneous and isotropic universe we must have 
Ea = 0 everywhere. Then equation o is violated, if pe 0. 

In fact, equation 0 holds in any spacetime if i?“ is interpreted as the electric field measured by an observer with 
a 4-velocity u“, pe as the charge density measured by the observer, and Da as the derivative operator associated with 
hab = dab + UaUb- In this sense, the classical equation ([2]) remains valid in general relativity, if it is hypothesized that 
the Maxwell equation m is correct. 
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Equation © can be integrated over a 3-volume V enclosed by a two-dimensional surface S' on E*. Then, by Stokes’s 
theorem, we get the Gauss’s law for electromagnetism 

£ J Ean^dS = Q , (8) 

where Q = fy pedV is the total electric charge contained in the volume V. If we apply the Gauss’s law in equation ([8]) 
to a volume in a homogeneous and isotropic universe, we get <5 = 0. Hence, Maxwell’s equations imply that inside any 
finite volume of a homogeneous and isotropic universe the total charge must be zero, conflicting with the assumption 
that the universe has a nonzero and uniformly distributed net charge. 

With the Gauss’s law in equation ([8]), a stronger conclusion can be obtained if the hypersurface Ej is compact (i.e., 
if the universe is closed). For a compact space Ej, we can use a compact surface S to divide it into two parts: Vi and 
V 2 , with V = Vi -I- V 2 . Both Vi and V 2 are three-dimensional compact manifolds with a boundary surface S. We can 
consider Vi as the interior of S, and V 2 as the exterior. Alternatively, we can also consider V 2 as the interior and Vi 
as the exterior (Fig. [T|). 

If the Gauss’s law in equation ([5]) is applied to the space Vi bounded by the surface S with an outward normal 
n“, as sketched in Fig. [U we get the total charge inside Vi: Qi = (l/47r) Jg EaU^'dS. If the Gauss’s law is applied 
to the space V 2 bounded by the surface S (then the is an inward normal), we get the total charge inside V 2 : 

Q 2 = — (l/47r) JgEan°‘dS = —Qi- Thus, Maxwell’s equations imply that the total charge in a closed universe must 
always he zero. This statement does not rely on the assumption of homogeneity and isotropy of the universe and 
uniform distribution of the charge. In other words. Maxwell’s equations cannot describe a closed universe with a 
nonzero net charge without any requirement on spacetime symmetry and charge distribution. 


III. MODIFICATION OF MAXWELL’S EQUATIONS 

In Sec. El we have shown that Maxwell’s equations cannot describe a universe with a nonzero total charge, which 
indicates that Maxwell’s equations have a fundamental flaw when they are applied to the universe as a whole. In this 
section and the rest part of the paper, we investigate how to modify Maxwell’s equations so that the inconsistency 
problem can be resolved. 

If we introduce an electromagnetic potential 4-vector A“ and define the electromagnetic field tensor Eab by 

Eab = ya^b — bAa , ( 9 ) 

the homogeneous Maxwell equation (U) is automatically satisfied and hence trivial. Since only the inhomogeneous 
Maxwell equation ([3]) leads to conflicting results when applied to a universe with a nonzero total charge, we need 
only modify equation (|3]). Since Fab = 0 in a homogeneous and isotropic universe, inclusion of additional terms 
proportional to Fab in the field equation does not help to address the problem. Then, the only possibility left is to 
add a term proportional to A“ to the field equation, if we want to keep the electromagnetic field equation linear in 
A“. 

The simplest way to modify the Maxwell equation m so that it can describe a universe with a nonzero total charge 
is to add a term —to the left-hand side, where ^ is constant. Then, equation ([3|) becomes 

= -47r. (10) 

When ^ is positive, this equation is equivalent to the Proca equation and can be interpreted 

as the rest mass of photons. Although ^ can just be a fundamental constant in nature, an effective term can 
naturally arise from spontaneous symmetry breaking by the Higgs mechanism if the electromagnetic field is coupled 
to a complex scalar field [H, . 

With the field equation the Type I Paradox described in Sec. El can be addressed. For instance, for an electric 
charge uniformly distributed in a homogeneous and isotropic universe, the solution to equation (HOI) is 

A“ = . (II) 

Here u°‘ = {d/dt)°‘ is the 4-velocity of an observer comoving with the expansion of the universe, and pe = Pe{t) is the 
charge density measured by him. The electric current density 4-vector is = PeU°'. 

By equation dm) we have 

A 'TT 

= — {paU‘^U^ + PeV^“) , 


( 12 ) 
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where Pe = dpe/dt. For a homogeneous and isotropic universe, is just the extrinsic curvature tensor of the 

hypersurface defined by t = const. Therefore we have and F^b = 0, which guarantees that the A°‘ is 

consistent with the cosmological principle and solves the electromagnetic field equation (fTOll . 

Contraction of any 4-velocity with equation (fTUl) leads to 

DaE- + = dTTPe , (13) 

where $ = —UaA°‘ is the scalar potential of electromagnetic fields. Integrating equation (IT^ over the volume Vi and 
V 2 in a compact space (Fig. [T]) respectively, we get 


Qi + Q 2 — f ^dV , (14) 

47r Jvi+V2 

which is not necessarily zero. Hence the Type II Paradox disappears too. 

The electromagnetic field equation (US can be derived from the action 

Sem = J (15) 

where the Lagrangian density Lem is defined by 

Lem = - ^FabF'^'^ - ^AaA^^ + J“ . ( 16 ) 

Here e is a fixed volume element, and g is the determinant of the spacetime metric in the coordinate system compatible 
with e. Variation of S'em with respect to Aa leads to the field equation 

Variation of S'em(>/“ = 0) with respect to g°'^ gives rise to the stress-energy tensor of electromagnetic fields [13 

The derived stress-energy tensor of electromagnetic fields described by the Lagrangian density in equation (ITBl) is 


where 


?EM,ab = ab + ^^^?EM,ab , (18) 

^°^TEM.ab = ^ (PacFb ^ J(19) 


is just the stress-energy tensor of electromagnetic fields in the ordinary Maxwell’s theory, and 

ab = ^ (^AaAb — -gabAcA’^ (20) 

is the additional stress-energy tensor associated with the ^AaA°' term in the Lagrangian. 

The trace of the electromagnetic stress-energy tensor defined by equations (fT51)-(E(Il) is Tem = —C^o^“/47r, which 
is nonzero unless is null. The divergence of the electromagnetic stress-energy tensor is calculated to be 

VaT|^ = -JaF^^ + , (21) 

47r 

where the field equation (nni has been applied. 

By the identity Vb^aF^-^ = 0, the field equation (ITOl) implies that 

Vaj“ = -fVaA^ . (22) 

47r 

When ^ ^ 0, the equation of charge conservation, Va = 0, is maintained if and only if 


Va4l“ = 0 , 


(23) 



which is just the Lorentz gauge condition. Substituting equation (l23ll into equation (|2T]) . we get the usual equation 
for the Lorentz force 


VaT^M = ■ (24) 

Equation indicates that when charge is conserved, the ^-term in the stress-energy tensor, i.e., the 
given by equation (1^ . has no effect on the dynamics of charged particles. However, the term ^^^TEM,ab affects the 
spacetime curvature through Einstein’s field equations. 

As an example showing that a photon mass term can arise from spontaneous symmetry breaking and the Higgs 
mechanism, let us consider a model with a Lagrangian density 

L = -jFabF-'^ - Vack - E(</)) . (25) 

Here, is a complex scalar field, and the asterisk denotes complex conjugate. The gauge covariant derivative operator 
TDa = Va + ieAai where e is the coupling parameter defining the interaction of the complex scalar field and the 
electromagnetic field. The potential function V{(/)) is taken to be of the form 

n<(') = |J-M"r</>+^(<(>*<(>)^ ( 26 ) 

where /r and A are positive and real-valued parameters. The constant term /r^/2A in V{(j)) is introduced to ensure 
that V{(f) = 0 at its minimum. 

The model described by equations ((^51) and (pS)) is the curved spacetime version of the model in the chapter 20.1 
of Ref. [ 3 ^ . The Lagrangian is obviously invariant under the gauge transformation defined by 

Aa —>■ Aa -t- Var(a;) , (27) 


and 




„ier{x),* 


(28) 


where r(a;) is any real-valued function. 

Since > 0, the potential E(</>) acquires a minimum at \(j)\ = (j)Q = Let us expand the Lagrangian about 

the vacuum state 4> = ^ 0 , and write (f) = 4 )q -\- [4>i{x) + i(j) 2 {x)]/\/2^ where and (j) 2 {x) are real-valued functions. 

The resulting Lagrangian is 

L = -\FabF^'^ - e^4>lAaA^ - - V^e^oA^^Vacj)2 + ... , (29) 

where “...” stands for terms cubic and quartic in Aa, ((> 1 , and 4 > 2 - 

For any given (j), we can always take the advantage of the freedom of gauge transformation to make (j) real. Hence, in 
this unitary gauge [11,133, we have 4)2 = 0. That is, the massless Goldstone bosons go away by a gauge transformation. 
Then we get a Lagrangian density which contains only physical particles 

L = -^FabF'^’^ - e^^AaA^ - ^VaC^iV^cj), - + ... , (30) 

where “...” stands for terms cubic and quartic in Aa and 0i. Clearly, the Lagrangian in equation (13011 describes a 
massive and real scalar field interacting with a massive electromagnetic field A“. The mass of the scalar field 4>i 
is m = v^/i. The mass of the vector field A“ (i.e., the mass of photons) is 

ruA = '/2e(f>o = ^[2 ep,jX^^^ . (31) 

Although spontaneous symmetry breaking and the Higgs mechanism naturally give rise to an electromagnetic field 
equation in the form of equation (HOD, they cannot be considered as the ultimate solution to the inconsistency problem 
in application of Maxwell’s equations to a universe with a net charge. When the symmetry is restored in high energy 
states of (j), the electromagnetic field becomes massless again and the inconsistency problem reappears. In the next 
section we will propose another type of modified electromagnetic field equations which, in our opinions, are better 
than the field equation and can be considered as a candidate for the ultimate solution for the inconsistency of 
Maxwell’s equations with a universe containing a net charge. 
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IV. EFFECTIVE ^-TERM FROM SPACETIME CURVATURE COUPLING 

Here, we propose a new type of modified electromagnetic field equations, where the electromagnetic potential vector 
is coupled to the Ricci tensor Rat of the background spacetime. We will show that the new field equations are 
better than the Proca equation and can also solve the paradoxes discussed in Sec. m 

Let us consider a Lagrangian density of electromagnetic fields with a term coupled to the Ricci tensor of the 
spacetime 


Lem = - '^RabA^A^ + Wa , (32) 

where k is a dimensionless constant of order unity and can be both positive and negative. The difference between 
this Lagrangian density and that in equation (fT6)l is that the {^/2)AaA°‘ term in equation ([T6|) is replaced by the 
{K/2)RabA°‘A^ term in equation (15^ . 

Variation with respect to Aa of the action with a Lagrangian density defined by equation (15^ leads to the electro¬ 
magnetic field equation 


- kR\A^ = -dTTj'’ , (33) 

which contains a term with the potential vector coupled to the Ricci tensor Rab- 

Let us check that equation (1551) can solve the Type I Paradox described in Sec. HT] For a net charge uniformly 
distributed in a homogeneous and isotropic universe, the cosmological principle requires that Fab = 0, = <I>(t)u“, 

and J“ = pe(t)u°‘, where = (i9/9t)“ is the 4-velocity of an observer comoving with the expansion of the universe. 
By the discussions in Sec. EH the potential vector = $(t)u“ automatically implies Fab = 0. For a homogeneous and 
isotropic universe, the Ricci curvature tensor Rab must have a form of Rab = a{t)uaUb+P{t)hab, where hab = gab+UaUb, 
a{t) and /3(t) are scalar functions of the cosmic time t. Therefore, we have R^aA°‘ = —since Uau‘^ = —1. 
Equation (1551) then has the following solution: $ = —{An/aiAiPe-, i.e., A°‘ = — (dTr/as:) 

Similar to the discussion in Sec. EH with the electromagnetic field equation (1551) . the Type II Paradox is also 
addressed. Contraction of any 4-velocity with equation (1551) leads to 

- KRabu'^A^ = 4^Pe . (34) 

Summation of the integrals over the volumes Vi and V 2 in a compact space (Fig. (H) leads to 

Qi + Q2 = -^ f RabU^A^dV , (35) 

dTT Jvi+V2 


which in general does not vanish. 

The electromagnetic field equation (1551) has the following features: 1. It does not introduce a new and small 
dimensional parameter, unlike the modified equation (1101) . It also does not rely on the existence of another matter 
field, e.g., a scalar Higgs field. The dimensionless k in equation (1551) is simply a numerical factor of order unity. 2. 
For a spacetime with Rab = 0, it becomes the standard Maxwell equation. This means that, unlike in the case of 
equation m, the modification to the Maxwell equation considered here does not affect the electrodynamics in a flat or 
Ricci-flat spacetime. In particular, in a flat or Ricci-flat spacetime photons remain massless and the gauge invariance 
is maintained. 3. For a spacetime with Rab = Apab, where A is the cosmological constant, the electromagnetic field 
equation (1551) reduces to equation (ITOl) with ^ = kA. 

Therefore, we consider the electromagnetic field equation (1331) more attractive and more fundamental than the 
Proca-type equation. If the universe contains a cosmological constant, by Einstein’s field equations we have 

Rab = Agab + Stt (^ab “ ]^Tga^ , (36) 

where Tab is the stress-energy tensor of matter, and T = g°‘^Tab- Then we get 

RabA^ = AAa + Stt (^ab " , (37) 

which indicates that the electromagnetic field equation contains a photon mass term with ^ = kA. In this sense, we 
can say that the electromagnetic field equation (1331) implies the Proca-type equation with a photon mass determined 
by the cosmological constant. 
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Application of Vf, to equation (l33)) leads to Va Jgff = 0^ where 


-'eff 


= J“ - 

47r 


Hence, the conservation of electric charge, VaJ°‘ = 0, implies that the gauge condition 

V“ {RabA^) = 0 


(38) 


(39) 


has to be satisfied. 

Derivations of the stress-energy tensor of the electromagnetic field from an action with the Lagrangian density in 
equation (|32l) . and the corresponding divergence, are presented in Appendixesand |B1 Here we only write down the 
results. The stress-energy tensor can be written as a sum of two terms as in equation m, with Tem, ab given by 
equation m, and 

^'^TEM.ab = ^ {V^^MaAb) - 2V'=V(,(A6)A,) + AA-R,(aAb) + 9ab [V,Vd(A^A'') - R.^A^A^] ] , (40) 

where parenthesis in the index of a tensor denote symmetrization of the tensor. The divergence of the stress-energy 
tensor is 

V“TEM.ab = ^ [Fbc - kR\A^) + (i?caA“)] . (41) 

47r 

Substituting the electromagnetic field equation (1^ and its divergence into equation (1^ . we get 

V“TEM.ab = -Fbar + AfeVa . (42) 

When the electric charge is conserved, we get the Lorentz force equation (EH). Hence, the electromagnetic force 
equation is not affected by the presence of a curvature-coupled term in the electromagnetic field equation. Physical 
meanings of this fact will be discussed in Sec. IVHl 


V. MOTIVATIONS FOR THE CURVATURE-COUPLED ELECTROMAGNETIC FIELD EQUATION 

Equation ([H is obtained from generalization of the Maxwell equation in a flat spacetime 

Q^F^^b ^ _^^jb ( 43 ) 

to a curved spacetime by direct application of the “minimal substitution rule”. That is, we replace in the equation 
the ordinary derivative operator da by the derivative operator Va associated with the metric tensor gab in a curved 
spacetime, and the F°'^ and J“ by their correspondences in a curved spacetime. In particular, the definition of Fab in 
a flat spacetime. Fab = daAb — dbAa, is replaced by the definition given by equation ()H- 

The “minimal substitution rule” provides a simple and convenient machine for converting an equation of physics in 
a flat spacetime to a curved spacetime. Although it can guarantee that the derived equation in a curved spacetime 
satisfies the principle of general covariance and returns to the correct equation in a flat spacetime, the “minimal 
substitution rule” cannot guarantee that the equation derived with it must be correct in physics in a curved spacetime. 
As a simple example, if the true field equation contains a term involving the spacetime curvature, the term is not 
present in a flat spacetime so cannot be recovered by the “minimal substitution rule”. A conformally invariant scalar 
field is just such an example (47j |. Since the term Rcjy/Q is missing in a flat spacetime, the equation of a conformally 
invariant scalar field cannot be derived from its version in a flat spacetime (i.e., the Klein-Gordon equation for a 
massless scalar field) by the “minimal substitution rule”. 

The “minimal substitution rule” has another profound problem. If a field equation in a flat spacetime contains 
second or higher order derivatives of a vector or a tensor, the corresponding field equation in a curved spacetime 
derived by the “minimal substitution rule” is not unique. Since the order of derivatives of a vector and a tensor 
matters in a curved spacetime, the “minimal substitution rule” can lead to different field equations in a curved 
spacetime even starting from the same equation in a flat spacetime, if the order of derivatives is arranged in different 
ways. The Maxwell equation (H51) contains second-order derivatives of the potential vector A“. So, even if the “minimal 
substitution rule” is insisted, different electromagnetic field equations can be obtained in a curved spacetime. 

To see the above point, let us write the Maxwell equation (l43l) in terms of A“ 

dad^A^ - dad^A^ = , 


(44) 
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which is equivalent to 


dad^A'^ - (1 + K)dad’’A‘^ + Kd^daA‘^ = 


for any k. Application of the “minimal substitution rule” to equation (|^5|) gives rise to 


Va(V“A^ - V^A“) - K (VaV^ - V^Va) A“ = -AttJ'’ , 


(45) 


(46) 


which is just equation (1551) by the identity = (VqV^ — V^Vq) A“. Therefore, since k can be any number (it 

can be any function in fact), the “minimal substitution rule” can give rise to an infinite number of field equations 
in a curved spacetime corresponding to the Maxwell equation in a flat spacetime. Although a requirement of gauge 
invariance picks up equation ([3]), an electromagnetic field equation in the form of (1551) cannot be excluded since gauge 
invariance is not a physical principle.^ We also remark that there is no observational or experimental evidence for 
gauge invariance of electromagnetic fields in a curved spacetime with a nonvanishing Ricci curvature. 

Although in equation (1551) it seems a little arbitrary in writing dad^A°‘ = (1 + K)dad^A‘^ — Kd^daA°‘ ^ there is a 
natural way to derive a field equation like that in (1331) . The value of k can even be determined. To do so, we first 
show that the Maxwell equation in a flat spacetime can be expressed in terms of a symmetric tensor instead of the 
anti-symmetric tensor. Let us define a symmetric tensor Hat by 

Hab = daAb + dbAa . (47) 


Then, by the identity 


Fab = Hab - 2dbAa , (48) 

the Maxwell equation (1451) can be written as 

daH^^ - d^H = , (49) 


where H = T]‘^^Hab = 2daA‘^. 

Now, applying the “minimal substitution rule” to equations (1571) and (1551) . we get the corresponding equations in 
a curved spacetime 


Hab — ^ aAb V bAa , 


and 


- V^’H = , 

where now H = g°'^Hab = 2VqA“. By the identity 

jdh 




equation (1551) is equivalent to 


VaF^^ + 2R\A‘^ = -Airr , 


(50) 

(51) 

(52) 

(53) 


which is just the equation (1551) with k = —2. 

The electromagnetic field equation m can be derived from an action with a Lagrangian density 


Lem = -^ {HabH^'^ - H^) + 4^Aa J“ , (54) 

by variation with respect to Aa- It can be checked that the above Lagrangian is identical to that in equation (1321) 
with K = —2, up to a boundary term which has no contribution to the action integral. So, the derived stress-energy 
tensor of the electromagnetic field is given by equations (1T8)) . (HH), and (150l) with tv = —2. 


^ According to Refs. [sSl and [S . gauge invariance is not a symmetry of nature. It generates nothing that is observable. While a global 
symmetry gives rise to a conserved current by Noether’s theorem, the local gauge symmetry does not. Gauge invariance only provides 
a principle for construction of a local theory for describing massless vector particles and a tool for the convenience of computations. 
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If we define another symmetric tensor 

Sab = Hab - Hgab , (55) 

the electromagnetic field equation (15111 can be recast in a neater form 

= -dTTj'' . (56) 

The Lagrangian density in equation ([54|) can be expressed with Qab as 

LeM = -1 (^0a60“'' - ^0') + Wa , (57) 

where 0 = g°'^Oab = —3i7. 

The above arguments support the proposal of equation (l33l) as a candidate for the electromagnetic field equation 
in a curved spacetime. The paradoxes arising from application of Maxwell’s equations to a universe as described in 
Sec.ini demands the inclusion of the curvature-coupled term in the electromagnetic field equation. The discussions 
in this section also explain why we consider to add a curvature-coupled term instead of a simpler RA^ to the 

field equation. 


VI. SOME CONSEQUENCES OF THE MODIFIED ELECTRODYNAMICS 


Both the two types of modified Maxwell’s equations discussed in previous sections have an unavoidable consequence: 
photons can have a nonzero rest mass. Equation (nni indicates an intrinsic photon mass, if ^ is interpreted as a 
fundamental constant. Equation (1331) indicates an effective photon mass as a function of the local spacetime curvature. 
Eor example, as mentioned in Sec. lIVl if there exists a cosmological constant as current observations indicated 
the kB!’^A°- term in equation (1551) gives rise to a photon mass = (zcA)^/^, corresponding to Aa 10^®cm (the 
size of the cosmic horizon). Here Aa = 27rAA = ‘li^h/mAC is the Compton wavelength of a massive photon. 

At present the most stringent and secure limit on the photon mass comes from the measurement of the magnetic 
field in the solar wind: toa < 2 x 10“®^g, or Aa > 2 x lO^^cm A much stronger but very speculative constraint 


comes from observations of the magnetic field on galactic scales, which is ttia < 10“°^g, or Aa > 3 x lO^^cm [^. In 
this section we investigate some consequences of a nonzero photon mass, or equivalently, of a positive Eor simplicity, 
we assume that vtia = is a constant. 


A. Propagation time delay of photons with different energy 

An immediate consequence of a photon mass is the appearance of dispersion in the velocity of photons: bluer photons 
travel faster, and the velocity approaches c as the frequency of photons approaches infinity. However, attempts in 
measuring or giving a limit on the photon mass by observing the difference in the arriving time of photons with 
different frequencies emitted by the same source are disfavored, since for any reasonable value of a photon mass the 
effect turns out to be extremely small 

For example, for a radio source at a cosmological distance D, the time delay of radio photons relative to a photon 
traveling with a speed c is 


At ( 




lOiOcm/ VlO^Hz 


-2 


where / = uj/2 tt is the frequency of a radio photon, and the cosmological distance D is defined by 

cdz 


Dx{z) = [ 
Jo 


H{z)il + zy 


(58) 


(59) 


Here z is the cosmic redshift, H = dlna/dt, where a = a{t) is the linear scale of the universe at time 

Taking Aa ^ 10^‘^cm, we get At < 10“®s for a radio source at a cosmological distance. This makes it almost 
impossible to detect the delay in the arriving time of photons due to a nonzero photon mass by observing the radio 
emissions of a gamma-ray burst or a quasar. 


® Tlie distance defined by equation ll59ll does not appear to be identical to any existing distance definition in cosmology, therefore we 
denote it by Dx- 
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B. Statistical mechanics of massive photons 

There is a profound discontinuity as the mass of photons approaches zero: massive photons have three degrees 
of freedom, but massless photons have only two. In this subsection we discuss the effect of a photon mass on the 
statistical mechanics of photons. To be specific, we assume that the background spacetime is a Minkowski spacetime. 

Massless photons obey the Bose-Einstein statistics and in the thermal equilibrium state have a spectrum of energy 
density given by the Planck radiation law 


gli a;3 
27r2c3 - 1 ■ 


(60) 


Here T is the temperature, k is the Boltzmann constant, and g is the statistical degeneracy factor. For massless 
photons we have g = 2. 

Massive photons are described by the statistics of relativistic Bose gases [5l,[53|. When kT ruAC^, the rest mass 
of photons can be ignored and the statistics approaches that of massless photons. The total number of photons are 
not conserved and the spectrum of energy density is given by equation O- The question is in the value of g. Naively, 
one may expect that g = 3 since massive photons have a longitudinal component (i.e., a spin component with s = 0) 
that massless photons do not have. However, the value of g depends on how the longitudinal photons interact with 
matter. If they interact with matter sufficiently weak or do not interact with matter at all, we still get g = 2 and the 
Stephan-Boltzmann law remains valid . 

Here we demonstrate that when photons have an energy hio ^ mA(? they interact with matter very weakly by 
considering the electromagnetic radiation generated by an electric dipole p oscillating with a frequency uj. Far from 
the dipole, the three-dimensional vector potential A and the scalar potential $ are respectively 


A(<,x) =-—, (61) 

r 

and 

$(t,x) = -—e-*‘^‘+*'='’p-n, ( 62 ) 

r 

where r is the distance from the dipole, n is the unit vector along the radial direction, and k = \JuP' — m\ is the 
magnitude of the three-dimensional wavevector. 

The magnetic field and the electric field derived from the A and $ are respectively 

B = ‘^e-^ku>+ikr (n X p) , (63) 

r 

and 

7 2 

E = -B X n -H ^ 

UJ r 

In the above expressions we have ignored terms of order and higher, since they do not contribute to the Poynting 
flux and the radiation power. 

The magnetic field B has only transverse components (i.e., B • n = 0). However, in addition to transverse compo¬ 
nents, the electric E has a longitudinal component E|| (i.e., E|| x n = 0) given by 

—iujt+ikr 

E|| = ^ ^ -(p-n)n. (65) 

From the expressions of B and E, we can calculate the time-averaged Poynting flux vector. When ^ = m\ ^ 0, a 
generalized Poynting flux vector is derived in Appendix IQ and given by equation (IC14|1 (see also [E^). Setting a = 1 
in equation (IC14I) , we get the generalized Poynting flux vector 

S = ^ [(5?(E) X K(B)) -I- m^(3?(<i>)3?(A))] = ^ sin^ 0 ^1 -I- cot^ 6 ^ n , (66) 


where () denotes time-average, and K(/) denotes the real part of a complex function /. As expected, S has only a 
component in the r-direction. 
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The power of the dipole radiation is given by the integration of S over a sphere of radius r, i.e., 


(67) 


As expected, the power is independent of the radius r as demanded by the law of energy conservation. 

Equation (1671) indicates that, in the radiation generated by an electric dipole, the power components in each photon 
spin state (s = 1,-1, and 0) have a ratio 


S=1 


: Pfs=-} '■ — 1 : 1 : 




( 68 ) 


Therefore, when w ^ ruA, the contribution of longitudinal photons to the total radiation is negligible. For example, 
if TTiA ^ 3.5 X 10“®2g (corresponding to Aa > lO^'^cm) and uj = lO^^Hz (corresponding to the peak frequency of 
the cosmic microwave background) , we get itiaIoo ^ 3 x 10“^® and < 10“^^. We cannot imagine any 

kind of electromagnetic radiation generated in lab or astrophysical conditions where longitudinal photons can have a 
detectable contribution. 

The relation in equation (1681) indicates that longitudinal and transverse photons are not statistically independent 
and do not have identical statistical weights. For a given statistical state, the number of longitudinal photons are 
always suppressed by a factor relative to the number of transverse photons. Because of the very weak 

interaction of longitudinal photons with matter, approaching to a statistical equilibrium state of longitudinal photons 
in a cavity is extremely slow, with a time-scale comparable to the age of the universe (sol - l^ . According to [6^, the 
practical impossibility for longitudinal photons to get equally partitioned with transverse photons in energy is due 
to the very large skin depth of longitudinal photons compared to that of transverse photons (with a ratio > w^/m^) 

[sllel. 

So, when kT ^ mAC^, we expect that for radiation containing both transverse and longitudinal photons, in the 
equilibrium state transverse photons have an energy spectrum determined by equation (ISOl) with g = 2, and the 
spectrum of longitudinal photons should be jointly determined by the spectrum of transverse photons and the power 
ratio relation given by equation (l68l) . Hence, for fio; ^ toac^, the energy spectrum of longitudinal photons should be 


u||(a;) 


m\d^ 

2n?uj^ 


u±(uj) 


m\c UJ 
2Tr'^h — 1 ’ 


(69) 


where Mj_(a;) is the energy spectrum of transverse photons (eg. 1501) . 

Finally, we argue that it is quite general that the contribution of m\ to the radiation field is ^ the 

dominant 1/r term as indicated in the dipole radiation case (^eas. l63l and l64p . In a Minkowski spacetime, from equation 
(HOD we can derive a generalized version of Ampere’s circuital law in the differential form 


V 


X B-— = 47r 

dt 



(70) 


Since dE/dt ^ —iwE and J = 0 outside the source, from equation (170)) we derive that 


E- (V X B-bm^A) . (71) 

UJ 

Since |V x B| ^ A:|B| and |B| = |V x A| ^ fc|A|, by equation (171]) we get |E| (from B) ~ (A:/a;)|B| ^ {k'^/uj)\A\, 
and |E| (from m^A) ^ (TO^/a;)|A|. Therefore we have 

|E| (from m\A) m\ 

|E| (from B) "" P ■ ^ ’ 

When kT <C ruAC^ we should have <C m\ and uj^ ~ Then, by equation (l68)) we get Ps=i ■ Ps=-i ■ Ps=o = 
1:1:1. Hence, when kT <C rriAC^, the statistics of massive photons become identical to that of a non-relativistic 
Bose gas with a degeneracy factor g = 3. As T —>■ 0, quantum effects become important and the Bose-Einstein 
condensation may occur [58l| . 


C. Constraint on the net charge in the universe and the charge difference between electrons and protons 

Assume that the universe has a mean net charge density 

Pe = xG^^'^Pm , 


(73) 
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where pm is the mean mass density, and the dimensionless parameter y is the ratio of charge to mass in the universe. 
The corresponding charge current density 4-vector is J“ = peu‘^, where u°‘ = is the 4-velocity vector of a fluid 

comoving with the expansion of the universe. For a homogeneous and isotropic universe, all quantities are functions 
of the cosmic time t. 

As discussed in Secs. m and uni in a homogeneous and isotropic universe we must have Fab = 0 everywhere. The 
solution to the electromagnetic equation (fTOjl is given by equation (fTTll . Then, by equations (fT^ - (l20l) . the stress-energy 
tensor of the electromagnetic field is 

T^m = jpI + ^5“") , (74) 


corresponding to a perfect fluid with mass density pA and pressure Pa defined by 

Pa = PAC^ = ^pI ■ (75) 

Substituting equation dzi into equation (HSl), we get PA = 2-kGc “^X^pIi^a, where Xa = ^ By pm = 

f4mPcrit(l + z)^, where z is the cosmic redshift, and pcrit = 1-9 x g cm“^ is the critical mass density of 

the universe, we get 

= ^ = A)’(! + *)», (76) 

Pcrit 4 \dH J 


Here the Hubble distance dn = c/FLq ~ 10^®ft.“^cm, and the Hubble constant Hq — 100/ikms“^Mpc~^. If the value 
of UJA at the present time (i.e., at z = 0) is denoted by we find loa = ^4^(1 + z'f. 

By the state equation (123, PA evolves with the cosmic redshift by pA ^ (1 + zXf’i i-C., pe oc (1 -I- z)^. Since 
Pra oc (1 -|- z)^, by equation (1731) the parameter y remains a constant as the universe evolves. 

Modern observations have confirmed that the matter content in the present universe is dominantly composed of 
about 70% dark energy (with P = —pc?), about 30% of cold dark matter and baryonic matter (with P <C p?), and 
a very small fraction of radiation (about 10“® in mass density with P = p?/5) [i^ Since the state equation 
Pa = PA? differs from all the matter contents just mentioned, comparison of pA with the mass density of various 
matter contents in the universe will provide a constraint on the value of y. 

For example, the requirement of pA < Pm at z = 0 leads to y < 2{iQ,m)~^^'^dH /Xa ~ /Xa- The requirement 

of PA < Pr (the mean mass density of radiation) at z = 0 leads to y < 2(Hr/3)^/^H“^d_f//AA ~ Q.QldH/XA- For 
comparison, y pz 10^^ and 10^® for an electron and a proton, respectively. 

A net charge in a universe can arise from several possibilities, including a difference in the fundamental charge of 
electrons and protons which together with neutrons form the fundamental blocks of ordinary matter, a nonvanishing 
charge of neutrons and neutrinos, and asymmetry of matter and anti-matter. Hence, a constraint on the parameter y 
can be converted to a constraint on the charge difference between electrons and protons, the charge of neutrons and 
neutrinos, and the asymmetry of matter and anti-matter 1^. [3^ . 

In the original paper of Lyttleton & Bondi [s^ , they (mistakenly) claimed that a small difference in the magnitude 
of electron and proton charges may provide a repulsive force for driving the expansion of the universe. Based on the 
then available observational data, they concluded that a fractional difference y ~ 2.2 x 10“^® would be sufficient to 
explain the observed expansion. However, as we have already demonstrated, in a homogeneous and isotropic universe 
Fab must vanish everywhere so there is no electromagnetic force on the cosmological scale. In fact, the stress-energy 
tensor of the electromagnetic field arising from the net charge causes an attractive force, instead of a repulsive force 
according to Einstein’s field equation. 

The approach of Lyttleton & Bondi can still be used to constrain the difference in the magnitude of electron and 
proton charges. An electron has a charge — |e|. Assuming that a proton has a charge (1 + y)\e\, where \y\ <C I. For 
a universe containing an equal number of electrons and protons, the total net charge density in the universe at any 
redshift z is given by pe = ypb\e\/mp, where rup is the proton mass, and pb is the average mass density of baryons in 
the universe. Comparison with equation (l73l) yields 


A = VXp 


fib 

flrr 


(77) 


where Xp = \e\IG^I'^mp = l.II x 10^® is the charge to mass ratio of a proton, and fib = Pb/Pait- Substituting equation 
(1771) into equation (1751) . we get 


WA = -^y Xpf^: 


1^] (1 + ^)® 


dn 


(78) 
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The requirement that uja{z = 0) = Ha < provides the first constraint on the value of y 


/Aa\ 

3i/2xpHfc [dn) 


1.3 X 10”^^ 



(79) 


where we have adopted flm = 0.3, Hb = 0.043, and h = 0.7 1^. 

A stronger constraint comes from the requirement that H j < H^.. Adopting = 4.2 x 10“^/i“^ (including the 
massless neutrinos and the cosmic microwave background [la[l3)i 


2Hy^ /Aa\ 

ydn) 


2.2 X 10"^® 



(80) 


An even stronger constraint on y comes from the condition of loa < Hr(l + z)'^ at z = 1100 (the redshift of cosmic 
recombination), where we have used the fact that the density of radiation oc (1 + z)'^. Since oja oc (1 + z)®, we get 


y < 


2ny^ 

3^/^Xp^b 



2.0 X 10"^® 



(81) 


The above constraint on the value of y is consistent with the result obtained in lab experiments if we take Aa dn 
10®®cm (sTjl . 

The estimate of the constraint on y presented above is based on the electromagnetic field equation m with a 
constant If ^ evolves with the cosmic time the above analyses must be accordingly modified. We also note that the 
above constraints are based simply on comparison of pA with the matter and radiation density in the universe. More 
complex analyses, e.g., comparison with the temperature fluctuations in the cosmic microwave background, may lead 
to stronger constraints on the charge asymmetry and the neutrino charge [H, 113 • 


VII. ON THE TESTABILITY OF THE FIELD EQUATION (f33|l 

In Sec. IIVI we have shown that when there is a cosmological constant A in the universe, the electromagnetic field 
equation (1331) contains a photon mass term with niA = (^A)^/®, corresponding to a Compton wavelength comparable 
to the Hubble distance. Given the very small value of a possible cosmological constant in the present universe [43,153, 
the effect of a so small photon mass may never be observable, unless the universe has a nonzero net charge ((Sec. IVip . 
However, in the epoch of inflation before the Big Bang when the expansion of the universe is presumably driven by 
a large effective cosmological constant, it can be imagined that a correspondingly large photon mass may affect the 
spectrum of the Gibbons-Hawking radiation [H] . 

Equation (1331) indicates that local curvature of the background spacetime can affect the electrodynamics. So, we 
can expect that, in a strongly curved spacetime, the effect of the curvature-coupled term may be measurable on scales 
comparable to the curvature radius of the spacetime. In this section, we estimate the order of magnitude of the 
curvature-coupled term, and find what kind of environment is favorable for detection of the effect. 

By equations (l33)l and (l38l) . the curvature term can be regarded as a pseudo-charge current density vector 

^ps = • (82) 

It is not a true electric charge current density, but its effect in generation of electromagnetic fields is equivalent to that 
of an electric charge current density. If an observer measures the electric field outside a charged ball and interprets 
the result with the standard Gauss’s law, he will find that the electric field appears as being generated by a total 
charge 


Qtot = Q + ^ [ RabU^A^dV . (83) 

47r Jy 

Here Q is the electric charge of the ball, is 4-velocity vector of the observer, and V is the volume of the ball. 

Similarly, if an observer measures the magnetic field outside a star and interprets the result with the standard 
Ampere’s circuital law, he will find that the magnetic field appears as being generated by a total current density 
inside the star given by 


jyt=r-^h‘^'RbcA\ 


(84) 
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where is the three-dimensional electric current density, and = 5 “*' -I- Therefore, even if the outside of an 
object is vacuum and Ricci-flat, the nonvanishing curvature-coupled term in the interior of the object can still affect 
the structure of the electromagnetic field outside the object. If the interior of the object is strongly curved due to the 
presence of a large mass density, the effect of the pseudo-charge current density arising from the curvature-coupled 
term may be testable in measurements of the outside electromagnetic field. 

To estimate the order of magnitude of the curvature-coupled term, let us ignore the cosmological constant and 
assume that the matter inside an object is described by a stress-energy tensor = pu°'u^, where p is the mass 
density of the matter. Then, equation dSIl) becomes 

RabA^ = 47rp{Aa-2^Ua) , (85) 


where d) = —UaA°‘ is the scalar electric potential. 

By Einstein’s field equations, the curvature radius of a mass system can be estimated by 


Zp 


( ^ 
\ PP 


- 1/2 


lO^^cm 


P 


Igcm 


-3 


- 1/2 


( 86 ) 


where Ip = 1.6 x 10 ^^cm is the Planck length, and pp = 5.2 x 10®^gcm ^ is the Planck mass density. Hence, inside 
an object, the order of magnitude of RabA^ can be estimated by where the curvature radius Vc is 

estimated by equation (|Mll . The first term in equation ([33|) can be estimated by |VaE“^| ^ \A^\/r\, where ta is the 
scale over which the electromagnetic field changes. Hence, we have 





[rj 


(87) 


Equation (l87l) indicates that the curvature term in the field equation (l33|) will lead to measurable effects if Tc < 
Objects clearly satisfying this criterion include neutron stars and the early universe. A typical neutron star has a 
radius r k, 10®cm and an overall mass density p r; 5 x lO^'^gcm”^, corresponding to a curvature radius ~ 10®cm 
by equation (1861) . Hence, we have ~ ta for an electromagnetic field with ta ^ t. Therefore, for a neutron 
star, we expect that the curvature-coupled term in the electromagnetic field equation can produce detectable effects 
on its inside and outside magnetic fields. In the very early time, the universe has a very high mass density and a 
very small curvature radius. The curvature-coupled term in the electromagnetic field equation should affect all the 
electromagnetic processes on scales of the curvature radius of the early universe. 

In the theory of electrodynamics defined by the Lagrangian in equation (1321) . the stress-energy tensor of electromag¬ 
netic fields contains an extra term in addition to the usual stress-energy tensor ^^^TEM,ab- The extra term ^-^^TEM.ah 
is defined by equation (l40)l . By equations (IB17I) and (1^ . the divergence of is 


V“(l)rEM.ab = FbaJps ” AfcVa Jp“ ■ (88) 

The ^^^TEM.ah only interacts with the pseudo-charge current. Hence, the extra stress-energy tensor ^^^TEM.ah represents 
a type of dark electromagnetic energy and momentum. It does not interact with electric charges and currents. 

In contrast, the stress-energy tensor ^^^TpM.ab leas. fT^ and IC4I) interacts with both the charge current J“ and the 
pseudo-charge current 7“^, since by equations (IB5I) and (1551) we have 

V“(°)rEM.ab = -Rba (+ Jps) ■ (89) 

In the total stress-energy tensor of electromagnetic fields, only the ^°^TEM,ab part is accessible to experiments based 
on electromagnetic methods. In a spacetime region that is strongly curved, the conservation of ^^^TpM.ab is affected 
by the presence of a pseudo-charge current, according to equation (1551) . 

Although l^lTEM.ab does not interact with electric charges, it affects the spacetime geometry through Einstein’s 
field equations. Note, ^-^^TkM.ab contains terms not involving the spacetime curvature (eq. 1401) . Even in a place 
where rA “C Cc (hence the curvature terms in eq. HHl are small), the magnitude of ^^^TEM,ab can be comparable to 
that of ^°^TEM,ab- In tbe early universe when the gravity of radiation drives the expansion of the universe, the dark 
electromagnetic energy may play a very important role. 


VIII. SUMMARY AND CONCLUSIONS 

Based on convincing arguments we have verified that Maxwell’s equations cause troubles and paradoxes as they 
are applied to the universe as a whole. Maxwell’s equations cannot describe the electrodynamics of a homogeneous 
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and isotropic universe with a uniformly distributed net charge, since in a homogeneous and isotropic universe the 
electromagnetic field must be vanishing everywhere (Type I Paradox). For a closed universe (i.e., a universe with 
a compact spatial section), Maxwell’s equations always fail if the total charge in the universe is nonzero, without 
any requirement on spacetime symmetry and charge distribution (Type II Paradox). We consider this issue as a 
fundamental flaw in Maxwell’s equations. 

We have investigated the possibilities for getting rid of the paradoxes by modifying Maxwell’s equations. Although 
an electromagnetic field equation containing a photon mass term (i.e., the Proca equation) can address the problem, 
the origin of the photon mass term itself need be explained. Therefore we have considered an electromagnetic field 
equation with a curvature-coupled term, i.e., the field equation (I33|) . The new field equation can also address the 
paradoxes and has the following advantages over the Proca equation: 

1. Unlike the Proca equation, the new field equation does not introduce a new dimensional parameter. In the Proca 
equation a photon mass is introduced as a fundamental constant. To be consistent with modern experiments 
and observations, the photon mass has to be incredibly small. In equation (I33L the /t is a dimensionless constant 
of order unity. 

2. In a flat or Ricci-flat spacetime, the new field equation returns to the Maxwell equation and gauge invariance is 
restored. In the Proca equation, the photon mass term is always there, and gauge invariance is always violated 
(though with introduction of a scalar field the gauge invariance can be recovered in the underlying theory). 

3. In a spacetime with a cosmological constant, the curvature-coupled term in the new field equation contains a 
photon mass term, where the photon mass is determined by the cosmological constant. In this sense, the new 
field equation implies the Proca equation and hence is more fundamental than the Proca equation. 

4. The new field equation has a strong motivation from the “minimal substitution rule”, as discussed in details 
in Sec. El It can naturally arise from the Maxwell equation when the Maxwell equation is generalized to a 
curved spacetime with the “minimal substitution rule”. The Proca equation has no similar motivation, unless 
the photon mass term is interpreted as arising from a more fundamental theory. 

5. The Proca equation may never be testable, if the presumed photon mass is incredibly small. The new field 
equation contains a curvature-coupled term, whose effect is determined by the strength of the Ricci tensor. In a 
spacetime region that is strongly curved, the effect can be important and measurable. So, the new field equation 
(1^ is testable. 

We have discussed some consequences of the modified electromagnetic field equations and the testability of the new 
field equation (1551) . Although in normal conditions the effect of a photon mass is hardly noticeable, a couple of places 
have been identified where the effect of the curvature-coupled term in the new field equation and the Lagrangian can 
be detectable. The places include the neighbor and interior of neutron stars, and the early universe. 

Due to the high mass density present inside a neutron star, a pseudo-charge current arising from the curvature- 
coupled term can affect the structure of the magnetic field of a neutron star. In the early time when the universe 
is small, the curvature-coupled term in the electromagnetic field equation affects all the electromagnetic processes 
on scales comparable to the curvature radius of the universe. The dark electromagnetic energy associated with the 
curvature-coupled term in the Lagrangian may also play an important role in driving the expansion of the universe. 

In summary, to address the inconsistency problem arising from application of the Maxwell equation to a universe 
with a net charge, we have proposed and investigated a new electromagnetic field equation (eq. 1551) . The new field 
equation is strongly motivated by consideration of generalization of the Maxwell equation in a flat spacetime to a 
curved spacetime, and is testable in a strongly curved spacetime region and in the early epoch of the universe. 
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Appendix A: Derivation of Equation (140 D 

By the equation 7.5.14 of [l^, we have 

2SRac = " g'’^VbVd5gac + 2g^‘^VbV^aSgc)d , 


(AI) 
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and 

2A‘^A‘=6Rac = - g^‘^A‘^A‘=VbVdSgac + 2g^'^A‘^A^WbWaSgcd ■ (A2) 

Let us evaluate each term on the right-hand side of equation (|A2p separately. For the first term, we have 

- g^^A-A^WaV,5gbd = -g^^^a {A^A^WJgbd) + 5“V,(AM^)V,(556d 

= -g“Va iA‘^A^VJgbd) + [Va{A<^A-)Sgbd] - g'^'^V,Va{A<^A<^)5gbd 
= V“[...]a - g^‘^V,Va{A-A-)Sgbd , (A3) 

where V“[...]a is not written out explicitly since it does not contribute to the integral of action so does not affect the 
derivation of stress-energy tensor. 

Similarly, for the second term on the right-hand side of equation (IA2I) . we have 

- g^^^A-A-S/bVdSgac = V“[...], - g^<^VdVh{A<^A-)5gac ■ (A4) 

For the third term, we have 

+ 2g'>'^A'^A-VbVJgcd = V“[...]a + 2g^‘^VaVb{A<^A-)5g,d ■ (AS) 

Hence we get 

2A-A-dRac = -g'>'^VcVa{A<^A-)5gbd - g^<^VdVM"A^)59ac + 2g'’^VaVfe(H“H^)dg,d , (A6) 

where = means “equal up to a term like V“[...]a”. Making use of the identity Sgab = —gacgbdSg‘^‘^, we get 

2A-A^dRac = [VeVd(A=H^)g„6 -b V‘'Vd(H„A) - 2V^Vb{AaA,)] 5g^^ . (A7) 

On the other hand, we have 

Rac5{A<^A-) = 2A^R,,Ab6g^^ , (A8) 

and 

RacA-A-dV^ = -^V^RcdA‘^A‘^gabSg‘^^ . (A9) 

Therefore, we have 

d {^gRacA^A^) = - V=V„(H,H,) + 2A<^R,,Ab + ^gab [VeVd(H^A‘') - R.dA^A'^] | 5g^^ . 

(AlO) 

Then, by equation (HU, the stress-energy tensor in equation (1401) is derived after symmetrization of the tensor index 
as required by the definition of Tab- 


Appendix B: Derivation of Equation (1411) 

For the simplicity of calculations, let us define 

nib = VeVd(A"A‘')g,b -f V^Va{AaAb) - 2V^V(„(Ab)A,) , 

and 

nf,) = 4A“i?,(„Ab) - RadA‘=A‘^gab ■ 

We have then 

TEM,ab = (°^rEM.ab + ^ (n^ + ni^.^) , 


(Bl) 


(B2) 


(B3) 
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and 


where 


V^TEM.ab = V“(°)TEM.a6 + ^V“ (h^ + 




EM,ab = -^FiiaVcF^ 
47r 


(B4) 


(B5) 


The divergence of 11^^^ is 

+ W‘^W‘^Wd{AaAb) - 2V“V=V(„(^)^,) 

= Vb (VeA'VrfA'' + A^VcVdA<^ + V.A'^VdA^ + A<^VcVdA^) + V“ (2V^A,VdA + + AbV‘^VdAa) 

- V“ (VMbVaA, + AbV^VaA, + V'AeVaA + - V“ (VM^Vb^e + A^V^V^A, + V=A,VbA, 

T^eV^Vb^a) 

= (2VbV,A"VdA‘^ + VbA^V.VdA'^ + A^VbV.VdA-^ + 2VbV,A‘‘VdA'^ + VbA‘'V,VdA‘= + A‘'VbV,VdA‘=) 

+ {2V^V‘^AaVdAb + 2V’^AaV^VdAb + V^AaV’^VdAb + AaV'^V’^VdAb + V^AbV’^VdAa + AbV^V^VdAa) 

- (y^V^^AbVaA, + V"AbV“VaA, + VMbV'VaA, + AbV“V"VaA, + V^V^A.VaAb + V"A,V“VaAb 
+V“A,V=VaAb + A,V“V"VaAb) - (V“V"A,VbA, + V"A,V“VfaA, + VMaV'VbA, + ^,V“V"VbA, 
+V“V"^eVbA, + V"AeV“VbAa + V“A,V"VbA, + A,V“V"VbA,) 

= Va X 0 + VbA'^ X @ + V^Ab X @ + X @ + [Rest] , (B6) 


where 


0 = 2 (VbVe - VcVb) = -2Rbd^^ @ = 0, 0 = (VaVe - VeVa)A“ = Rcddl^ (B7) 

0 = 2(VbVe - VeVb)dl“ + (VcV“ - V“Ve) Ab = -2R,^a “A'" + i?bdc "A^ , (B8) 


and 


[Rest] = A^Vb^c^dA‘‘ + A^Vb^c^dA^ + A,V<^V<^VdAb + AbV’^V'^VdAa - AbV“V=V„A, - A^V^V^V^Ab 
-A,V“V=VbA, - A,V“V^VbA„ 

= A" [Vfa (V,VaA“ + VaVcA“) - (VcVa + VaVj VfcA^] + Ab (VaV, - V,Va) V"A“ 

+A“ (VaV, - VcVa) V"Ab . (B9) 

By the equation 3.2.12 of [l^, we have VaVbWcd — VbVaWcd = Rabc + R-abd^^ce for any tensor ujab, from which 
we get book'd - = -Rabe’^^^'d + Rabd^'^’^e- Hence, we have 

(VcVa + VaVj V6A“ = 2V,VaVbA“ + R,,, "VeA“ + ReeVfaA" , 


(BIO) 


where we have used the identities R„ 


_ _ TD a _ _ TD 

— — Jtrp* 


We also have VcVaA“ + VoVcA“ = 2VcVaA“ + RcdA‘^. Then, in the expression for [Rest], 

A" [Vb (VcVaA“ + VaVcA“) - (V,Va + Vc) VfcA^] = A" [2Vc (VbVaA“ - VaVbA“) + Vb (RcdA'^) 
-R,,b "VeA“ - ReeVbA'^] = A^ [-2Ve (RbdA'") + Vb {RcdA‘^) - R^eb "VeA“ - R.eVbA^] , 

where we have used the identity VbVcVaA“ = VcVbVaA“. 

For the other two terms in [Rest], we have 


Ab (VaVe - VeVj V"A“ = Ab = -AbRae (V^A“ - V“A") = 0 , 


and 


A“ (VaV, - V,Va) V=Ab = A“ {-RaeV^Ab + R„,b "V=Ae) . 

Therefore, we have 

[Rest] = A= [-2V, {RbdA'^) + Vb (RcdA'") - R,,, ^VeA“ - R.eVbA^] + A“ {-RaeV^Ab + R„,b "VM, 

= A" [-2Ve (RbdA'") + Vb (RcdA'") + Reabe(V^A“ + VM") - R,e (VbA'^ + V^Ab)] . 


(Bll) 

(B12) 

(B13) 

(B14) 
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Putting the above results together, we get 

V“nW = _2Ve [A^RmA^) + A^Vb (RcdA^^) - A'^RdcVbA^ , (B15) 

where we have used the identity -Rbcda + Rbdca + Rdcba = 0 (since R^^^cf = 0)- 
( 2 ) 

The divergence of Il^f, is 

= 2Ve {A^RdbA^) + 2V^ {A^RdcAb) - Vb (RcdA^A^) . (B16) 


Therefore, we have 

V“ = 2V“ (A^RdcAb) - 2A'^RdcVbA'^ = -2Fb'^R,dA<^ + 2AV" [R^dA'^) . (B17) 

Substituting equations (IBs!) and (IB 171) into equation (jB4l) . we get equation (HTl) . 

Note, when we derived the equation (IB17|) we did not make use of the electromagnetic field equation. When 
Rab = 0, we have 11^^^ = 0 and = 0. 


Appendix C: Derivation of the Poynting Flux Vector 

In this appendix we derive the Poynting flux vector corresponding to the electromagnetic field equation (nnj. 

By equations ([5]) and (|6]), the Fab can be expressed in terms of Ea and Ba as 

Fab = —2.Es^aUb] — CabcdB'^u’^ , (Cl) 

where = —1 and EaU°‘ = BaU°‘ = Q. We can derive 

EabF^'^ = -2 {EaE^ - BaB^) , (C2) 

and 

FacEb^ = —EaEb — BaBb + EaE'^UaUb + {gab + UaUb)B(.B‘^ + 2u(^a^b)cdE'^B^ , (C3) 


where iabc = -^abcdu'^- 
By equation (IT^ . we get 

47r^'’^rEM,ab = {EcE"^ + BcB'^) ^UaUb + ~ EaEb — BaBb + 2U(^ah)cdE‘^B^ . (C4) 

To express the law of energy conservation in terms of Ea and Ba, we assume that the spacetime has a timelike 
Killing vector parallel to the i.e., = au“, where a is the lapse function. By the Killing equation we get 

0 = V(o4) = aV(an&) + U(oVb)a , (C5) 

from which we can derive 

aVau“ + = 0 (C6) 

and 

u“VaMf) = Vf, Ina — lira . (C7) 

Since lA’ub = —1 and aUb = 0, equation (IC7I) leads to 

M“Vaa = 0 , u“VaWf) = Vb Ina . (C8) 


Then, by equation (|C6p we have 


Vau“ = 0. 


(C9) 
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An energy-momentum flux vector is defined by 

Ja = —■ 
By equations (IT^ . (IC4I) . and (l20l) . we get 

a f 1 ■ 


ja ^ 


4^ { i \e,E^ + Bc-B" + C ($" + icA")] u“ + + e$i“) | , 


(CIO) 


(Cll) 


where we have written A“ in terms of a scalar potential $ and a spatial vector potential A“ (i.e., A“ = -|- A“, 

= 0 ). 

With the expression in equation (ICllI) . the energy-momentum flux vector is decomposed into an energy density 
component parallel to u“, and a spatial momentum component perpendicular to it“: 


= PEMU^ + T 


where 


and 


Pem = 


a 

Stt - 


E^E^ -h -b ^ -b AcA' 


Att V 


-b ) . 


(C12) 

(C13) 

(CM) 


In the first term is just the usual Poynting flux vector, and the second term is an additional momentum flux 
due to the ^-term in the electromagnetic equation. We may call the generalized Poynting flux vector. 
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